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Abstract. We calculate the time- and energy-dependent 
opacity of high-energy 7-rays attenuated by pair- 
production interactions with accretion-flare photons that 
are scattered by gas and dust surrounding the nuclei of 
active galaxies. We show that the temporal behavior of 
the high-energy opacity cutoff can be used in conjunction 
with the time history of the accretion flare to determine 
the location of the 7-ray emission site, as well as the col- 
umn density and scale height of the material surrounding 
the central engine. Reverberation mapping using this tech- 
nique is now possible for nearby BL Lac objects such as 
Mrk 421, and will be particularly valuable when the 10 
GeV - 100 GeV 7-ray window is opened, because here the 
attenuation of distant blazar radiation by the intergalactic 
infrared background radiation is negligible. 

Key words: Galaxies: nuclei - galaxies: active ~ gamma 
rays: galaxies - black hole physics 



1. Introduction 

Markarian 421, a BL Lac object with redshift z = 0.031, 
is the first extragalactic object detected at TeV energies 
with high significance (Punch et al. 1992). Searches for 
other AGNs emitting VHE 7-rays have been conducted 
using the air Cerenkov telescopes of the Whipple Obser- 
vatory (Schubnell et al. 1994). These searches have nat- 
urally focused on those extragalactic 7-ray sources al- 
ready detected at photon energies i? > 100 MeV with the 
EGRET instrument on the Compton Gamma Ray Obser- 
vatory (Fichtel et al. 1994; Dermer & Gehrels 1995), but 
have so far yielded negative results. Two explanations have 
been advanced to explain these observations. Because Mrk 
421 is so close by cosmological standards, its detection, 
coupled with the failure to detect more distant blazars, 
implies that TeV 7-rays from the more distant objects 
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are attenuated by the diffuse intergalactic infrared radia- 
tion field (Stecker, de Jager, & Salamon 1992; de Jager, 
Stecker, & Salamon 1994). By contrast, Dermer & Schlick- 
eiser (1994) point out that BL Lac is a truly lineless object 
(see, e.g., Kikuchi & Mikami 1987), and as such has only 
very tenuous emission line clouds near the central engine 
compared to the stong emission lines found in quasars 
(and some BL Lac objects). Thus the rescattered central 
source radiation is very dilute in Mrk 421, permitting the 
TeV 7-rays to escape without attenuation. 

In 1994 May, Mrk 421 was observed to flare in TeV 
7-rays during a period of ~ 2 days by a factor of 10 
over its flux level in the period 1994 Jan. -Apr., or by 
a factor of 5 above its quiescent flux level measured in 
1992 and 1993 (Kerrick et al. 1995; Macomb et al. 1995). 
Observations with ASCA made 2 days after the onset of 
the TeV 7-ray flare also showed an enhanced level of 2- 
10 keV X-rays compared with earlier iJAOS'^ T" measure- 
ments (Takahashi et al. 1994). These observations suggest 
a method to diagnose the structure and density of the gas 
and dust surrounding the central source of Mrk 421 and 
other 7-ray blazars. By monitoring the time dependence 
of the flaring emission from the central source and the pair 
attenuation cutoff of the high energy 7-rays due to photon- 
photon pair production, we can extract information about 
the column density, number density, and location of the 
site of high-energy 7-ray production. This is a new ex- 
ample of a reverberation mapping technique which uses 
correlated temporal signatures to monitor conditions of 
a physical system, analogous to methods where variabil- 
ity of emission line intensities are monitored in response 
to the temporal behavior of photoionizing radiation from 
the central source (cf. Peterson et al. 1991). 

In Section 2, we describe the system and model it un- 
der the assumption that the central source radiation can 
be approximated by a point source. Asymptotic forms for 
the total photon number density and opacity at times long 
after the peak of the flare are derived in Section 3, and 
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numerical simulations are presented in Section 4. We sum- 
marize in Section 5. 



2. Time-dependence of photon opacity from scat- 
tered radiation 

Attenuation of MeV-TeV 7-rays by photons with energies 

in the infrared through X-ray range can occur in the cen- 
tral regions of active galactic nuclei if both the intensity of 
the lower energy radiation and the scattering depth of the 
surrounding medium are sufficiently great, provided that 
the spatial extent of the scattering medium is not too large 
(e.g., Blandford 1993; Dermer & Schlickeiser 1994; Bland- 
ford & Levinson 1995). The presence of strong emission 
lines from quasars implies the existence of surrounding 
clouds of gas which can rescatter the central-source ra- 
diation (Sikora, Begelman, & Rees 1994). The structure 
and column density of the surrounding gas near the su- 
permassive black holes powering quasars is, however, not 
well known. Information about the location of the 7-ray 
emission site and the spatial extent and density of gas 
and dust surrounding the supermassive black hole can be 
extracted using the method we now describe. 

We consider the attenuation of 7-rays emitted outward 
along the z-axis at height Zi above a central soft-photon 
source (e. g., an accretion disk) which exhibits a flare de- 
scribed, for simplicity, by a Gaussian in time: 



Nph{e, f2; t) = Nmax{e, ^) e 



(1) 



In eq. (1), Nph,{e,fl;t)d€dSl is the differential number of 
photons emitted per second at time t with dimension- 
less energy e = hv/uieC^ in the range e and e + de that 
are directed into the solid angle element dCl in the di- 
rection ri = {0,4>). The full- width half maximum dura- 
tion of the flare is given by 2^/^(ln2)cr = 1.96c7, and the 
photon emissivity at the peak of the flare is denoted by 
N„iax{^,^)- Some of these soft photons are scattered by 
surrounding clouds of density ne(r) into the path of the 
7-rays. The scattered soft photons are directed at an an- 
gle 6* = cos^^ /i* with respect to the z-axis, as shown 
in Figure 1. We assume that the gas is moving nonrela- 
tivistically and can be treated as a stationary scattering 
medium. We also assume that the gas is optically thin 
to Thomson scattering so that only single scatterings are 
important. 

The opacity of a photon with dimensionless energy ei 
due to photon-photon pair production can be written as 



00 27V 1 00 

Tj^{ei,ti,Zi) = J dz J dtp* J dn*{l - fi*) J 

Zi -1 



de- 



ei(l-(<*) 



Fig. 1. Geometry of the system. Soft photons are emitted at 

time t from an accretion disk which is approximated by a point 
source at the origin. The photons are Thomson-scattered at 
location r by electrons in surrounding dust and gas clouds. 
The scattered radiation intercepts high-energy photons which 
are beamed outward along the z-axis, and contribute to the 
opacity of the high-energy gamma rays through 7 — 7 pair 
production attenuation. 



(Gould & Schreder 1967; for corrections see Brown, 
Mikaelian, & Gould 1973). Here ti denotes the emission 
time of the high-energy photon, nph{e, jj.* , (p*; z,t) repre- 
sents the spectral density of scattered soft photons per 
unit solid angle at height z and time t, and the pair pro- 
duction cross section is given by 



•[(3-/3^)ln([i|) -2/5(2-^2)] 



(3) 



where 



eei(l-M*) 



(4) 



■aj^{ei,e,fj,*) Hphie, ij,* ,(/)*■, z,ti H ^). (2) 



is the velocity of the electron/positron in the center-of- 
momentum frame. 

The soft photon density per unit solid angle at location 
r in the direction Q = {0, </>) defined by the location of 
the scattering material relative to the position of the soft 
photon source is given by 
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nph{e,n;r,t) 



S{n - n). (5) 



We approximate the Thomson scattering event by 
isotropic scattering with cross section ctt and neglect re- 
coil, which is a good approximation when e <^1. The pho- 
ton production rate by Thomson scattering at location r 
and time t is therefore 



(It 



^{v)Nph{e„il;t--^ , (6) 



where the subscript "s" is used to denote the scattered 
quantities. 

Because the particles at each point in the cloud are 

assumed to scatter isotropically, the spectral density of 
scattered photons at height z along the z-axis and at time 
t can be obtained by integrating the scattered spectrum 
over all space. We find that 



1 oo 

/Th (r) • / z — Z' — V — X 
dr -^Nph ie,fl;ti + 



-1 



<5(M*-r) 

(10) 



In order to perform the /i* -integration by evaluating 
the (5-function it is necessary to change the order of inte- 
grations. The limits for the remaining integrals then follow 
from the threshold condition 



eei(l-r) >2, 



(11) 



and we find 



oo oo 1 oo 

T~f^{ei,ti,Zi) = ^ J dz J de J djj, J ^^^^^(1-/2*)- 

Zi A. —k Tri 



Z — Zi — T — X 

•f7T,^(e, ei,/!*) Nph ( e,fJ;ii H ) , (12) 



where 



d^jdfijdr 
-1 



ph 



r + X 



(7) 



where we have dropped the subscript denoting scattered 
quantities, and the overbar on Q in the expression for Nph- 
It is then a simple matter to obtain the angle-dependence 
of the photon spectral density through the expression 



i oo 

nph{t,^i*,(t>*;z,t) = ^ J dn j dr 



"-e(r) 



-1 



r + X 



s{r-r)s{i^*-r{r,fi,x]). (s) 



Here x = {r"^ + z'^ — 2r2/i)^/^ is the distance between the 
scattering electron and the z-axis at height z, the angle Jl* 
is given by 



(9) 



and the angle (f>* = tt + (see Fig. 1). Substituting equa- 
tion (8) into equation (2), assuming azimuthal symmetry 
of the soft-photon emission spectrum and distribution of 
scattering material, gives the photon opacity 



rr...=zV^H ''^,-_lf^ ) (13) 
if /i^ ^ fc^ , and 



A; = 1- 



eei 



(14) 



If /i = fc wc have Vmin = ^ (for this case /c > is re- 
quired), while for /X = — A; condition (11) can not be satis- 
fied {rmin oo). 

To illustrate the importance of the effect of photon 
attenuation by reverberating soft photons, we assume that 
the disk radiates isotropically and that its radiation during 
the flare is described by a power law in photon energy 
described by the expression 



Nph{e,n;t) 



Kt 



Air 



(15) 



with cutoff energies e' and e". The power-law spectra could 
be produced in the inner hot regions of an accretion disk 
by nonthermal electrons which radiate synchrotron pho- 
tons at frequencies above the self-absorption frequency. A 
more detailed model of the flaring emission from an accre- 
tion disk would require time-dependent solutions of the ac- 
cretion disk equations (see, e.g.. White & Lightman 1990), 
which we neglect for simplicity. We also assume for pur- 
poses of illustration that the surrounding clouds are ho- 
mogeneously distributed with uniform density ne(r) = Uq- 
The constant K is related to the maximum luminosity 
of the disk by 
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where L46 is the maximum disk luminosity in units of 10^^ 
ergs s~^ and mgC^ is the electron rest mass energy in ergs. 



nph{e,iJ.*, (j)*;z,t) {2Tr) ^nph{e; z,t)6{fi* - Hq) , (21) 



3. Asymptotic behavior of photon density and 

opacity we find 

If we consider the effects of reverberation at times much 

later than the flare maximum, considerable simplifications r-y^(ei, ti) = ^""^'^"^'^ (ctj — Zi) ■ 

of the exact expressions (7) and (12) for the scattered 

photon density and opacity, respectively, arc possible. The 

derived asymptotic expressions provide useful checks on °? 

the numerical results presented in the next section. • / -r— / de e~" cr^^(ei, e, /xJJ) 

Under the assumptions made at the end of the previous /. * 

section, we may write the total scattered photon density 
as 



ei(l-»'n) 



2noKaTcr /ei\"-i f dz 



00 r+z 



(t-s±i)" V27f V2y J (cti + z-Zi) 

, . noaTKe~°' /" , /" , e 2^ 

nph{e;z,t) = — dr dx . (17) 

8ncz J J xr 

\r-z\ 



At times a and t^-, equation (16) reduces to 



a °° °° (^-xq)^^ where 
nph{e;z,t) J — dr J dx ^ 



(l-/^Sr y rf/?/?(l-/32)"-2a^^(/3) (22) 



-00 



•e(r + 0-a;o)e(xo-|r-^|), (18) \v otherwise 



where x^ = ct — r and G denotes the Heaviside function 
such that 0(y) = 1 if 2/ > and 6(t/) = otherwise. Thus, 

we have 



and 



noarK e"" ^ ^ f dr 1 Here we may assume « 0, (3rnax « 1 for the model 

n'ph(^', z,tj w Sttcz ^ ^''^^ / — _ ^ -^ve consider, because the pair production cross section (3) 

vanishes for /? — !■ and for /3 ^ 1. 

The 2-integration in eq. (22) can easily be done ana- 
no arK e~'^ \/2 7rcr^ ^ ( ct + z\ lytically, and the second integration is only dependent on 



Attcz ct \ct—z) the spectral index a. Thus, we define 



noCTTKe " \/27ra^ , , 



1 



(^*)' ' I{a)= fd00{l-0^r-^a^^i0) (24) 

Note that this last expression does not explicitly depend ^ 
on z. 

Inserting this expression into eq. (2), we set t = ti + -^(1-5) ~ 0.3 ar 

, noting that approximation (19) is valid for ct^ Zi. 
The disk photon density at the location of the scattering The function I{a) is only weakly dependent on a. The 
electrons decreases as ^ and the scattered photon density, photon opacity for times late after the flare is given by 

again, decreases as -3-. Therefore, the main contribution 

of scattered photons will come from that angle 9* with 

respect to the 2:-axis for which = r. This yields / 2^/^ no i^'or c 1 1 rt \ 

-'^-^^^'^ = ^17^ 

= (20) (25) 

CI, i- z which is independent of Zi because we neglected zi com- 

Using the scattered photon density pared to cti in equation (25). 
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4. Numerical Results 

We now show the results of numerically performing the 

integrals (7) and (12), and compare these results with ap- 
proximations (19) and (25), respectively, derived in Sec- 
tion 3. We assume that the medium surrounding the cen- 
tral engine has a Thomson scattering optical depth of 0.02, 
extends out to a scale height of 0.1 pc, and is uniformly 
distributed in this region. This yields an electron den- 
sity no = 10"'' cm^'^. Furthermore, we describe the disk 
photon spectrum during the flare by a power law with 
spectral index a = 1.5 extending from dimensionless pho- 
ton energy = 10~^ (corresponding to 5 • 10""^ eV) to 
e" = 10~^ (5 • 10^ eV), and assume its maximum lumi- 
nosity Lmax = 10*^ ergs s~^. 



the approximation because the latter does not take into 
account the lack of scattering material at r > 0.1 pc. 



Fig. 3. Time dependence of the pair-production opacity 
for photons emitted at different heights Zi given by the labels 
on the curves. Other parameters are the same as in Fig. 2. The 
dotted curve shows approximation (25) for r^-y. 



Fig. 2. Time dependence of the scattered photon number den- 
sity and pair-produetion opacity t~,~/. The photon number den- 
sity is computed at z = 10^^ pc, and r-y-y is for a 500 GeV 
photon (ei = 10®) interacting with lower-energy photons pro- 
duced in an accretion-disk flare described by a Gaussian in 
time with a full-width half-maximum duration of one week. 
The approximations valid at late time are shown by the dot- 
ted curves. 

The time-dependent eS'cct of a one-week flare de- 
scribed by equation (15) on the pair-production opacity 
for 500 GeV photons (ei = 10*^) emitted at Zi = lO'^ 
pc, together with the total scattered-photon number den- 
sity at the injection height, is shown in Fig. 2. We see 
that the photon opacity reaches maximum values of or- 
der unity and, due to the flnite light travel time, attains 
its maximum at times later than at t = Zi/c which, for 
the example considered here, is ~ 10^ s. The numerical 
results show good agreement with the approximations for 
the asymptotic behavior at late times, which indicates a 
decrease of the opacity as 1/t. If ct > 0.1 pc {t ^ 10^ s), 
of course, the numerical results will decrease faster than 



In Fig. 3 we show the dependence of the photon opac- 
ity on the injection height Zj, with all other parameters 
the same as given in the case shown in Fig. 2. From vari- 
ability time scale arguments, Kerrick et al. (1995) derived 
an upper limit of 0.03 pc on the size of the emitting re- 
gion of the TeV radiation for the case of Mrk 421. This 
is in accord with the values of Zi used in Fig. 3 if the in- 
jection height and inferred size scale are comparable. As 
the injection height increases, the maximum photon opac- 
ity decreases and attains its greatest value at increasingly 
later times. It is interesting to note that, as we already 
saw from equation (25), the photon opacity late after the 
flare maximum does not depend on the injection height. 
Note also that the delay between t = Zi/c and the time of 
maximum opacity depends very weakly on the injection 
height. In the case considered here, the delay At w 10^ s, 
corresponding to the duration of the flare. 

The energy dependence of the opacity T-y-^ is shown 
in Fig. 4 where, again, the dotted curves represent the 
analytic approximation (25). For these calculations, we 
assumed an injection height of Zi = 10^^ pc. From this 
flgure, we see that the form of the time-dependent photon 
opacity is csentially the same at all gamma-ray energies. 
The relative amplitude of the opacity curves is described 

by 

T^^(ei,ii)(xer' , (26) 
as was shown for the asymptotic behavior (eq. [25]). 
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Fig. 5 shows the time-dependent modification of a pri- 
mary 7-ray spectrum by the effects of pair-production 
opacity effects from accretion disk-flare photons. Here we 
let the primary gamma-ray spectrum be described by a 
power law with spectral index a = 2.0, and choose an 
emission height Zi = 10~^ pc. The results are plotted for 
different times with respect to the time of the maximum 
of the accretion-disk flare. As can be seen, 7-rays with 
energies ^ 500 GeV will be almost completely absorbed 
after ~ 10^ s following the flare maximum and will remain 
strongly attenuated long after the flare. If we assume a 
greater injection height, this effect will occur correspond- 
ingly later and the influence of reverberating photons will 
be restricted to higher energy 7-rays, as can be seen from 
Fig. 4. 



Fig. 4. Time dependence of the pair-production opacity Tj^ 
for photons emitted at height Zi = 10^"' pc with dimonsionloss 
energies given by the labels on the curves. Other parameters 
are the same as in Fig. 2. The dotted curves show approxima- 
tion (25) for T-y-y. 



Fig. 5. Time dependence of the observed photon spectrum 

due to photon-photon attenuation by soft photons from an 
accretion-disk flare which are scattered into the high-energy 
photon beam. Here we assume a power-law photon spectrum 
with photon spectral index a = 2. The high-energy photons 
are assumed to be emitted at height Zi = 10~* pc. Other par 
ramctcrs arc the same as in Fig. 2. The curves are labeled by 
time after flare ma^ximum. 



5. Discussion and Summary 

We have considered a highly simplified geometry to il- 
lustrate the potential of gamma-ray astronomy to map 
the structure and density of clouds of dust and gas sur- 
rounding the central nucleus of a BL Lac object or quasar. 
Photon-photon attenuation intrinsic to the source will oc- 
cur only if the surrounding medium has a sufSciently large 
column density or if the gamma-ray emission site is s\if- 
ficiently near the central nucleus. As pointed out in the 
Introduction, the scattering emission line clouds are rel- 
atively tenuous in Mrk 421 compared to quasars. It is 
therefore likely that the number density of the scattering 
medium in Mrk 421 is significantly less than 10^ cm~^, 
the value assumed in our model calculations. As the pair 
production optical depth scales roughly linearly with the 
number density of the scattering material, this is consis- 
tent with the observational result by Kerrick ct al. (1995), 
who did not find any temporal variation of the spectral 
shape during the TeV 7-ray flare in 1994 May. 

This observation implies a limit on the density of 
the scattering material surrounding the central engine of 
Mrk 421 if the TeV 7-ray production site is within 0.03 
pc (see Section 4). Using Eq. (26) and the results shown 
in Fig. 3, and assuming that the pair production optical 
depth up to several TeV (ei 10^) always remained less 
than 0.1, we find that 



no ;i lO'' LJe cm"^ w 3 L^^ cm-=* (27) 

for a Gaussian flare with a FWHM duration of 1 week. 

Optical and TeV gamma-ray campaigns to monitor Mrk 
421 can therefore provide constraints on the density of the 
surrounding medium of Mrk 421. During the May 1994 
outburst of Mrk 421, however, only variations in the keV 
and TeV ranges were observed, whereas no significant vari- 
ations were seen in near-simultaneous mm, IR, UV, and 
GeV gamma-ray observations (Macomb et al. 1995). The 
keV X-rays are most effective in attenuating 0.1- 1 GeV 
gamma rays. We could apply our results to this regime, 
noting that the X-ray flare luminosity was at the level of 



-1 cm-3 



3L; 



-1 cm-3 
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~ 2 • 10^^ ergs s ^. Although the poor statistics of GeV 

observations might not support a constraint as strong as 
0.1, equation (27) should yield a sensible estimate 
for this case. This implies a density of no ^ 10^ cm~^ for a 
one week X-ray flare, corresponding to a Thomson scatter- 
ing optical depth of tt ~ 0.002 of the surrounding scatter- 
ing medium if it exends out to a scale height of 0.1 pc. Ac- 
counting for probable inhomogeneitics (c. g., a decrease of 
the cloud density with distance from the center) together 
with the relatively high location of the 7-ray emission site, 
this limit on the Thomson scattering optical depth should 
be raised somewhat in order to yield a more conservative 
estimate. It should be kept in mind, however, that the X- 
ray flare in this event is probably beamed emission from 
the jet, in which case these estimates do not apply. 

As shown in Fig. 5, the time-dependent truncation of 
the gamma-ray spectra following an accretion-disk flare 
gives useful information about the location of the gamma- 
ray emission region and the mean density of the surround- 
ing gas and dust. By comparing the temporal behavior of 
the flare and the attenuated gamma-ray spectrum in de- 
tail, further information about the distribution of scatter- 
ing gas can be obtained. 

We can predict idealized time-dependent behavior of 7- 
ray emission in response to variations in the lower-energy 
radiation for specific models. In the models of Dermer & 
Schlickeiser (1993) and Sikora et al. (1994), gamma-ray 
production in blazars is caused by Compton scattering 
of accretion-disk photons or rcscattcrcd radiation. Flar- 
ing accretion-disk radiation in the infrared to X-ray range 
will produce flaring behavior at 100 MeV to GeV ener- 
gies due to the time-variable soft photon energy density 
in the relativistically outflowing fluid frame. If there are 
extremely energetic electrons in the fluid frame, then a 
flare can also occur in the TcV range. The shape of the 
unabsorbed IC-spectrum has been calculated for isotropic 
relativistic electrons in the outflowing fluid frame (Dermer 
& Schlickeiser 1993). If the accretion-disk flaring behavior 
in the infrared and optical range shows the same behav- 
ior as the UV/X-ray flare, then the 7-ray spectrum will 
be absorbed according to the time and energy dependence 
illustrated in Fig. 5. 

It is important to note that it is the infrared and opti- 
cal photons that arc most effective in attenuating TeV 
gamma rays (e.g., Stecker et al. 1992). To apply our 
method to TeV measurements, therefore, correlated in- 
frared/optical and high-energy gamma-ray observations 
are necessary (note that no space-based detectors may be 
necessary in this case). However, the diffuse intergalactic 
infrared background radiation produced by the superpo- 
sition of reprocessed emission from stellar nucleosynthe- 
sis renders most quasars with redshifts ^ 0.2 invisible, 
due to strong 7 — 7 attenuation over this pathlength (see 
also Stecker & de Jager 1993; MacMinn & Primack 1994). 
For applications to a larger number of objects, it is there- 
fore desirable to devise 7-ray detection systems in the 10 



GeV - several hundred GeV regime, where the effects of 

the diffuse infrared background radiation arc small. The 
identification of attenuation intrinsic to the source will 
be necessary, in any case, to implement methods to infer 
the Hubble constant from high-energy gamma-ray obser- 
vations (Salamon, Stecker, & De Jager (1994)). 

Finally, we note that the presence of scattered radia- 
tion in the source will confuse any attempt to determine 
the level of the intergalactic infrared background radiation 
using high energy 7-ray observations of blazars (Stecker 
et al. 1992; Dermer & Schlickeiser 1994). By monitoring 
temporal variations at low and high-photon energies, our 
method provides a way to subtract out the level of atten- 
uation intrinsic to the source, and determine whether the 
diffuse intergalactic infrared background or rescattered 
radiation is more important for attenuating high-energy 
photons in a given source. 

Acknowledgements. We thank Dr. R. Schlickeiser for comments 

and discussions. Useful comments by Dr. John Mattox and the 
referee, Professor O. C. do Jager, arc also acknowledged. Par- 
tial support for the work of C. D. was provided under NASA 
grant DPR S-30931-F. M. Bottcher acknoledges financial sup- 
port by the Deutsche Forschungsgemeinschaft. 

References 

Blandford, R. D. 1993, in Compton Gamma-Ray Observatory, 

ed. M. Friedlander, N. Gehrels, and D. J. Macomb (New 

York: AIP), p. 533 
Blandford, R. D., and Levinson, A. 1995, ApJ, in press 
Brown, R. W., Mikaelian, K. O., & Gould, R. J. 1973, Ap 

Letters, 14, 203 
De Jager, O. C., Stecker, F. W., and Salamon, M. H. 1994, 

Nature, 369, 294 
Dermer, C. D., and Schlickeiser, R., 1993, ApJ, 416, 458 
Dermer, C. D., and Schlickeiser, R., 1994, ApJS, 90, 945 
Dermer, C. D., and Gehrels, N., 1995, ApJ, in press 
Fichtel, C. E., et al. 1994, ApJS, 94, 551 
Gould, R. J., and Schreder, G. P. 1967, Phys. Rev., 155, 1404 
Kerrick, A. D., et al. 1995, ApJ Letters, 438, L59 
Kikuchi, S., and Mikami, Y. 1987, PASJ, 39, 237 
MacMinn, D., & Primack, J. 1994, preprint 
Macomb, D. J. et al. 1995, in preparation 
Peterson, B. M., et al. 1991, ApJ, 369, 119 
Punch, M., et al., 1992, Nature, 358, 477 (1992) 
Salamon, M. H., Stecker, F. W., & De Jager, O. C. 1994, ApJ, 

423, LI 

Schubnell, M., et al. 1993, in The Second Compton Sympo- 
sium, ed. C. E. Fichtel, N. Gehrels, and J. P. Norris, (New 
York: AIP), p. 597 

Sikora, M., Begelman, M. C, and Rees, M. J. 1994, ApJ, 421, 
153 

Stecker, F. W., and de Jager, O. C. 1993, ApJ, 415, L71 
Stecker, F. W., de Jager, O. C, and Salamon, M. H. 1992, 

ApJ, 390, L49 

Takahashi, T., et al. 1994, lAU Chcular No. 5993, 1994 May 
23 

White, T. R., & Lightman, A. P. 1990, ApJ, 352, 495 



8 



M. Bottcher, C. D. Dermer: Reverberation Mapping of AGN 



This article was processed by the author using Springer- Verlag 
M^X A&A style file L-AA version 3. 



n I I r 



scattered photon number density [e " cm ^] 



(approximation) 




scattered photon number density [e " cm ^] 

at z = 10"^ pc 



Photon number at location of the disk [6*10^^ e " s ^] 



T (approximation) 




(numerical) 



J L 



1 



8 9 10 11 



t[10 s] 



10 



10 







z. = 10 pc 



Disk photon number 
atz = 

[6*10'" e'^s"'] 



z. = 10"" pc 




z, = 10"^ pc 



10 



-1 



10" 
t[s] 



10' 



